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I. INTRODUCTION 



Transitions among different phases have been observed as temperature increases or de- 
creases, in wide variety of physical systems from elementary particles to condensed matters. 
They have attracted great interest both experimentally and theoretically, and profound 
physics content has been revealed particularly in infinite systems. Despite an isolated sys- 
tem, statistical properties of atomic nuclei have been investigated in terms of the so-called 
nuclear temperature providing us with a ground to study phase transitions in nuclei. 
However, it is not straightforward to identify phase transitions in finite systems like nuclei, 
because the quantum fluctuations associated with finiteness of the system obscure the tran- 
sitions 0, 0| so that even definition of 'phases' should depend on theoretical models to a 
considerable extent. 

It will be natural, from theoretical points of view, to define phases in nuclei via the mean- 
field picture, since 'phases' in physical systems are more or less a semi-classical concept. 
Indeed, we discuss spherical/deformed and normal-fiuid/superfiuid phases in nuclei based 
on the mean-field approximation. Phase transitions are often connected to the spontaneous 
symmetry breaking, leading to violation of conservation laws. In the mean-field picture of 
nuclei, the rotational symmetry is broken in the deformed phase, and the particle number 
conservation is violated in the superfiuid phase. However, the symmetry breaking in nuclei is 
a hypothetical effect due to the mean-field approximation. In practice, the conservation laws 
are restored via a sort of quantum fluctuations arising from correlations beyond the mean 
field. In studying statistical properties of finite systems, it will be significant to investigate 
effects of the symmetry restoration. 

As the most general mean-field theory, we shall primarily consider the Hartree-Fock- 
Bogoliubov (HFB) approximation. To restore the symmetries, projection with respect to 
the quantum numbers should be applied j4|. The particle number projection as well as the 
angular momentum projection methods have been developed for zero-temperature cases; i.e. 
for the ground-state wave function. Bi, and also to. the quantum-n^mbe.-constrainad 
HFB (CHFB) solutions along the yrast line 0, B 1^- For finite-temperature problems, 
projections with respect to the discrete symmetries such as the parity and the number-parity 
were explored more than two decades ago 01 • However, there is a certain complication in 
extending them to the continuous symmetries, arising from the non-commutability between 
the projection operator and the Boltzmann-Gibbs operator in the mean-field approximation. 
A way to overcome these problems has been found in Ref. jll| for the particle number 
projection within the Bardeen-Cooper-Schrieffer (BCS) approximation. General projection 
formalism in the variation-before- projection (VBP) scheme is developed in Ref. and is 
extensively applied in the framework of the static-path approximation. 

In this paper we reformulate the quantum number projection at finite temperature by 
employing the thermo field dynamics (TFD) 0, Q|. The TFD has been shown to be a 
powerful tool to handle the thermal fluctuations within the mean-field theories. We 
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shall show that this is also true for the quantum-number-projected HFB theory at finite 
temperature. While the resultant formulae are equivalent to those in Ref. the present 
formalism is advantageous in the following respects. In the TFD, the thermal expectation 
value of an observable is expressed in terms of a vacuum expectation value in an enlarged 
TFD space. Therefore explicit formulae are derived in a straightforward manner, keeping 
complete parallel to the zero-temperature formalism We demonstrate it for the simul- 
taneous projection both of the particle number and the angular momentum. This leads to 
another advantage that variational parameters can be identified in an obvious manner and 
handled easily. While application of the present method to the VBP calculations is straight- 
forward, further approximation is desired in order to carry out the variation-after-projection 
(VAP) calculations. We shall also discuss an approximation of the entropy, so as for the 
Peierls inequality hold which guarantees the variational principle for free energy. It 

is expected that this approximation scheme will make the VAP calculations practical. 

In Sec. ini application of the thermo field dynamics (TFD) to the HFB theory at finite 
temperature is presented. In Sec. IIIIl we formulate the TFD version of the quantum-number 
projection in the HFB at finite temperature. As well as general arguments, explicit formulae 
for the particle number and angular momentum projection are presented by taking specific 
basis-sets. Formulae in the BCS approximation are also given. In Sec. IIV| we propose an 
approximation scheme for the entropy, which keeps the Peierls inequality. Finally, the paper 
is summarized in Sec. IVl 



II. FINITE-TEMPERATURE HFB THEORY AND THERMO FIELD DYNAM- 
ICS 

A. HFB equation at finite temperature 

In the variational derivation of the Hartree-Fock-Bogoliubov (HFB) equation at finite 
temperature 10, a set of variational parameters {Uk^,Vki_i,Ul^,V^^} is introduced 
through the generalized Bogoliubov transformation (GBT), which relates the original single- 
particle (s.p.) operators {cfc,c^} to those in the quasiparticle (q.p.) picture {afj,,ajj}, 



The transformation matrix W obeys the unitarity relation W'^W = WW^ = 1. We express 
the original s.p. state as \i) = c||vac.), with the vacuum |vac.) satisfying Cj|vac.) = 0. The 
q.p. vacuum |0) oc Il^au\va.c.) is annihilated by |0) = 0. Note that (0|0) = 1. 

An additional set of variational parameters {E^j^} comes into the problem through the trial 
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statistical operator 



Wo 



; Ho = Y^ E^ala^ , (2) 



Tr(e-^o/^) 

where we put the Boltzmann constant k-Q = 1. We denote the thermal average of an operator 
O by (O) = Tr(ii)oO), where the trace is taken over the grand canonical ensemble. 
We consider the Hamiltonian comprised of up to two-body interactions, 

-f^ = ^ eic\ci + ^ Vijkic\c]ciCk, (3) 

i ijkl 

where two kinds of nucleons are not discriminated to avoid unnecessary complication. The 
hermiticity of the Hamiltonian implies Ei = e* and Vijki = —Vjiki = —Vijik = vl^^. In the 
constrained HFB (CHFB) approximation at finite temperature, the grand potential is given 
by 

= {H') -TS; H' = H-\Z- KN - uj,,J, , (4) 

where Ap (An) and c^rot are Lagrange multipliers, which are interpreted as proton (neu- 
tron) chemical potential and rotational frequency of the system, respectively. The proton- 
(neutron-) number operator Z (N) and the x-component of angular momentum operator Jx 
are expressed in terms of the s.p. operators as 

^ = Y1 ^i^*' ^ = Y1 ^^i^'^ ■^'^ = {'i'\Jx\3)c\cj, (5) 

iep ien «i(all) 

where the sum Ylii^p (^ien) extends over the proton (neutron) states, and ^jj^^u) over 
both proton and neutron states. In the CHFB the conservations of proton number, neutron 
number and angular momentum are taken into account in their averages. 



{Z) = Z, {N) = N, {.l) = ^JiJ+l), (6) 

which indicate constraints. In calculations without some of the constraints, the associated 
Lagrange multipliers are set to be zero in Eq. 

The auxiliary Hamiltonian in Eq. is rewritten in terms of the q.p. operators by the 
GBT dH), 

H' = Uq + Hii + H20 + H22 + H^i + H^Q. (7) 

For later convenience, explicit forms in the above expression are presented in Appendix A. 
Ensemble averages of the bilinear forms in the q.p. operators are given by 

{a^Uy) = {a^al) = , (8b) 
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where is the occupation number of the /i-th q.p. level. Thus, the variation of the 
parameter 6E^ is converted to 5/^ = — /^(l — f^)5E^/T. The approximate entropy S in 
Eq. (jH) is expressed as 

S = -Tr(t&olnt&o) = - J][^ln/^ + (1 - /^)ln(l - f,)] . (9) 

The s.p. density p and the pair tensor k, are defined by 

p,, = (c]q) = - f)V'' + UfU%, , (10a) 

= (c,q) = - f )U'' + t//\/t]^,. . (lob) 

The ensemble average of H' becomes 

{H') = Trs.p. Up + ^Tp + =Uo + Y.{H,,),,f, + 2 Y,{H,,),,,J,U , (H) 

where Trg.p. stands for the trace over the s.p. space, 

^ij = {£i - K)Sij - UJrot{i\Jx\j) (12) 

with r = p (n) for proton (neutron), while the HF potential matrix F and the pair potential 
matrix A are defined by 

Tij = {T^)ij = y^^VikjiPik , (13a) 

kl 

Aij = -Aji = ^"^Vijkitiki ■ (13b) 

kl 

The variational principle 5^1 = 6{{H') — TS) = yields 

(Hlf)^^ = [U^{^ + T)U -V^i^ + Vyv + U^AV -V^A*U]^ 



I pv 



p 

dS 



(14a) 



J pu 



(^2o)m^ = [u\^ + T)v* -v^^ + ryu* + U^AU* -V^A*V*]^ 

= iH2o)p. + 6 5^(if3i)M.pp/p = . (14b) 
p 

These equations are summarized in the form of the HFB coupled equation at finite tem- 
perature. 

In the CHFB case this equation should be solved together with the constraints in Eq. ©, 
by which the Lagrange multipliers Ap, An and t^rot, as well as self-consistent q.p. energies 
E^, are determined as functions of quantum numbers and temperature. 
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B. Bogoliubov transformation extended by TFD 



In the thermo field dynamics (TFD) |l3, 14|, the original q.p. operator space {a^, aj^} is 
enlarged by including newly introduced tilded operators; {a^, a^, aj^, aj^}. Correspondingly, 
the q.p. vacuum in the enlarged space is defined by |0) ® |0), where |0) is the tilded vacuum 
satisfying q;^|0) = 0. The TFD vacuum |0) ® |0) is hereafter denoted by |0) for the sake 
of simplicity, i.e. a^\0) = a^\0) = 0. In order to handle ensemble averages in a thermal 
equilibrium at temperature T, the TFD vacuum |0) is transformed to the temperature- 
dependent vacuum by the following unitary transformation. 



|Or) = e-^|0), G = -G^ = J2^ 



(16) 



where 



sin ■}}. 



(17) 



This unitary transformation relates the operator set {a^, a^, , aj^} to new set of the 
temperature-dependent operators {Pfi, Pfj., Pl, Pl}, 



a. 



a, 



Gn -G 



e^^e-^ = g^P^ - g^, 



(18) 



so that the temperature-dependent vacuum should fulfill [3^\Qt) = P^j\^t) = 0. The 
temperature-dependent quasiparticles created by or will hereafter be called TFD 
quasiparticles. The thermal average of an observable O is then expressed by the vacuum 
expectation value at |0r), 

(O) = Tr(wod) = (0t|O|0t). (19) 

Now we unify two unitary transformations and ()18p to compose an extended form of 
the GBT Q, 



f I3\ 

P' 

VP' I 



Wi 




(20) 



where 



kfj, 



Ug V*g 
-V*g Ug 



kfi 



kfi 



Vg U*g 
-U*g Vg 



kfi 



(21) 



Because of the unitarity relation WW'^ = W'^W = 1, the matrix inverse to W is given by 



W~ 




(22) 
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III. QUANTUM NUMBER PROJECTION IN FINITE-TEMPERATURE HFB 
THEORY 



A. Projection operators 

1. General form 

Since quantum numbers are usually associated with a certain group structure of the 
system, projection with respect to the quantum numbers is also introduced in connection to 
the group. In general, a projection operator P^u) is defined as a subring basis corresponding 
to an irreducible representation (irrep.) f?*^"^ of a group Q jiol. liot. 

^■' = ^^^^E^£'(--')-. (23) 

^ x&g 

where x stands for an element of ^, dim(f)'^"^) is the dimension of the representation matrix 
Q^°'\ and g the order of Q. The property of group representation yields 

pl:^p^' = , (24) 

justifying a convenient bracket representation of the projection operator, Pj°^ = |a/i)(az/|. 
The projection on the subspace corresponding to the irrep. without referencing /x, is 
obtained by 

^« = E^M?=El"'")("'"l- (25) 

It is obvious that Pa is idempotent, i.e. P^ = Pa- 

When Q is not simple and is decomposed into a direct product of an invariant subgroup 
and the complementary quotient group as Qi ® its irrep. g^"^ is also a product such as 
® f?2"^\ where f?-"''* is an irrep. of Qi {i = 1,2). Denoting the projection operator on 
f?!"'"* by Pat, the projection operator on g^"'^ is Pa = Pa^Pa2- Therefore, projection operators 
of simple groups are essential. 

In the following we assume ^ to be a compact simple group or a product of such groups, 
whose element x is represented by an appropriate unitary operator Q = e~*'^. Here S = S{Q) 
is a hermitian operator, which belongs to a Lie algebra associated with Q and is dependent 
on a set of parameters 0. The projection operator Pa is represented in the integral form 



Pa = J dQCa{e)Q, (26) 
where Co;(Q) is an appropriate function of and is derived from Eq. 
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2. Angular momentum projection operator 



For the angular momentum projection, the relevant group is SU{2), which is denoted by 
SU{2)j in this paper. The group element is the rotation operator R, whose parameters are 
represented by In practice, we consider two alternative parameterizations. One is the 
Euler angles $ = (a,/3, 7). The rotation is also represented by an angle u {0 < u < 2tt) 
around an axis indicated by a unit vector n. Using the representation = sin 9 cos 0, Uy = 
sin^sin0,n^ = cos 6', we arrive at the other parameterization $ = {00,6, (p). Corresponding 
to these parameterizations, the rotation operator is expressed in two ways as 

= Q~^f^-^v = Q-'i^i'n-xJx+nyJy+n^J^) (27) 

where the angular momentum operators are defined in terms of the s.p. operators as in 
Eq. ©, 

Jk= ^{AJk\j)c\cj {k = x,y,z). (28) 

ii(all) 

The angles {u,9,(f)) are related to the Euler angles by 

a; (3 a + 7 (3 / . a + it + a - 

cos — = COS — COS , tan u = tan — / sm , = . 29 

2 2 2 2/ 2 ^ 2 ^ ^ 

The order of the group is determined by the volume of the parameter space, 

g= da sin f3df3 / d-f = A sin^ -rftu / sin^d^ / rf0 = 167^^ (30) 
Jo Jo Jo Jo 2 Jq Jq 

where the range of 7 is taken to be < 7 < 4-7? for the double-valued representation 

corresponding to half-odd-integer spin. The irrep. of SU{2)j is given by the Wigner D- 



f unction 
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Di^{^) ^ {JM\RmJK) = e— ^4fi.(/5)e-^^^ , (31) 

which is a unitary matrix of dimension (2 J+1). Then Eq. (j23|) yields the following projection 
operator. 



PiiK = = ^TTF da fAnPdp d^ D';,j,{^)R{^). (32) 

-LOTT Jq Jq Jq 

For integer spin J, the integral J^^ d^y may be replaced by 2 ^7. 
We express the character of the representation matrix Dj^j^i^) as 

xAu) - D^^($) - . (33) 

M=-J ^ ' ' 



221. If 



The real function Xj('^) can be represented in many different ways as listed in Ref. 
the magnitude of angular momentum J is subject to the projection while its z-component 
in the laboratory frame M is not referenced, the projection operator is given by 

2 7-1-1 /"^'^ in 1"^ r'^'^ 
Pj = P]^ J2 HiM = ^^ sm'-du sinede d<PxAu;)m). (34) 



M=-J 
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3. Particle number projection operators 

The relevant group to the particle number projection is the gauge group f/(l), whose 
group elements are parameterized by a single variable ip for the particle number operator 
N; exp{—iipN). For neutrons, the group is denoted by U{l)iy, and the projection operator, 
which projects out states having the exact neutron number A^, is given by, 

PN = Ph = ^ r e-^^"(^-^)rf¥^n. (35) 

^TT Jo 

Likewise, the proton number projection is implemented by the operator, 

= Pt = i_ (36) 

and the relevant group to Pz is denoted by U{l)z- The product operator PzPn is employed 
for the U{l)z x f/(l)jv projection, by which the canonical trace is calculable from the grand- 
canonical trace. 



4- SU{2)j X U{l)z X U{l)iy projection operator 

The projector of simultaneous projection of angular momentum J, proton number Z and 
neutron number is P(j^z,n) = PjPzPn, which satisfies PfjzN) ~ P{J,z,n)- This P(j,z,n) 
can be expressed in the form of Eq. (pUj) . with 6 = ($, ipp, ipn)- The operator Q stands for 

Q = P($) e-^'^p^e"*'^"^ = e-'^, (37) 

where 

S = (ppZ + {priN + un ■ J; n ■ J = UxJx + ^yJy + n^Jz ■ (38) 

Obviously P{j^z,n) commutes with the total nuclear Hamiltonian H, which conserves angular 
momentum, nucleon number and charge. However, note that [P(^j^z,n), Hq] ^ [l^ . 



5. Parity and number-parity projection operators 



The space reflection forms the discrete group isomorphic to 5*2 . The projection operator 
with respect to the space parity can be expressed in the form similar to Eq. ()26|) . with the 
integral converted to a discrete sum. 



1 + ^exp {iiT^clck) 



(39) 
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Here ^ = +1 (.^ = —1) for the projection of positive (negative) parity state, and denotes 
the sum extending only over the s.p. states of negative parity. As far as the parity is not 
mixed in the q.p. state fi, the projection operator is also written as 



1 + ^exp (ivr^aj^a^) 



(40) 



as given in Ref. 

The number-parity is relevant to another 5*2 group, which is a discrete subgroup oiU{l)z 
or of U{1)n- Respective to protons and neutrons, the number-parity projection is carried 
out by the operator 



1 + 1] exp [in ^ clck) 

k 



(41) 



where rj = +1 [rj = —1) for the projection of the state of even (odd) number-parity. Apart 
from the denominator, this is obtained by restricting ipj^ in P/v (or ipp in Pz) only to and vr. 
Since the GBT does not mix the number-parity, Pg is equivalently represented in terms 
of the q.p. operators [l^ . 



1 + 7] exp (ivr^aj^a^) 



(42) 



B. Extended form of linear transformation 



We now consider the projected statistics. The projected ensemble average of an operator 
O is newly defined by 

Tr{Pe-P^op) Tr(e-/^^oP) 
where the statistical operator with projection is introduced as 

Wp = . (44) 

Tr(e-/3^op) 

For an operator commutable with P, (0)p is calculated in the TFD by 

(0)p = MwoOP) ^ (0t|OP|0t) .45. 
Tr(u;oP) (Ot|P|Ot) ' 

Substituting the projector by the form of Eq. ()2(i|l . we obtain 

/6)p = im{OT\dQ\OT)de 
/C(e)(Or|Q|OT)rfe ■ 
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Thus, in course of projection calculations, we meet with the TFD vacuum expectation values 
of following types: 

1 

/5l/32/33/?4 

where /3 represents any of the TFD q.p. operator of four types, or The 

operator Q is given in Eqs. (|37l38p for the SU{2)j x U{l)z x f/(l)jv projection. 



(0^ 



Q|Or), 



(47) 



We here define the transformation matrix of the s.p. operator with respect to the unitary 
transformation Q, assuming that Q conserves the particle number. 



kk' cy- 

This linear transformation is extended to the TFD space by 

/cfc\ ( Qkk' \ 



(48) 



Q 



Ck 

'^k 



( Ck' \ 



Ck' 
Ck' 

V4/ 



(49) 
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\ 1 / 

Applying the TFD-extension of the GBT in Eqs. ()20I21I22|) . we convert the above transfor- 
mation to the one among the TFD q.p. operators. 



Q 



k,k' V 




(Qkk' 

1 
Ql^, 
\ 

( 

Pi 






1/ 



k'v 



(50) 



In the last expression, there appear four types of the matrices K, L, M and N, which play 
central roles in the projection calculation 0,01 • present formalism, these matrices, 

whose dimension is also doubled, are given by 



T* 



M, 



N. 



giU^QU + V^Q*V)g + giU^QV* + V^Q*U*)g 
giy^'QU + WQ*V)g g{V''QV* + W'Q*U*)g + g^ 

giV^'QU + U''Q*V)g g{VQV* + WQ*U*)g - gg 
giWQU + V^Q*V)g - gg g{WQV* + V^Q*U*)g 



where we have employed abbreviated notations, e.g., 

m^QU + V^Q*V)g + g\, = UUl.Qkk'Uwu + V:^Qkk'*Vk',)g^. + 



(51a) 
(51b) 

(52) 



k,k' 
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C. Wick's theorem and reduction formulae 



Corresponding to Eq. 



the operator S in Q = e has the form 
S = Skk'clck' ■ 



(53) 



k,k' 



According to the general prescription of projection calculation 16|, we represent S in terms 
of the TFD q.p. operators as 



(54) 



/x,i'(all) 



/^,!^(all) 



where the summation ^^^^^^^ extends over both the tilded and non-tilded q.p. states. On 
the other hand, the TFD-extension of the GBT derives an alternative form as 



Sii Si2 

S21 S22 



22 )lJ-f 



( I3u\ 

PI 
\PI) 



+ 



'I, ] \Ol2)^,u 



Pi 
~Pl 



+ h.c. 



(55) 



where ^n, 5*12, 5*21 and 6*22 are calculated as 

iS22)f_iu = 



{gV'^SUg),, {gV'^SV*g),, 

{gV''SV*-g)^, (gV^'SUg),, 
{gWSV*g)^, (gWSUg)^, 

{gWSV*g)^, (gWSUg)^, 
{gV''^SV*g)^, {gV'^SUg)^, 

Comparing ()54p with the last expression in ()55|) . we obtain 

5^°) = Tr^22 , S^^~^ = Su - S^2 , S'-^^ = 2^12 = 2Sl 

By applying the general formalism of projection the TFD vacuum expectation value of 
the operator Q is given by 



'2\ 



(56a) 
(56b) 
(56c) 

(57) 



(OtIQIOt) = (detir)^/2exp 



^(0) + ^TrS^i) 



(detir)^/^exp 



(68) 



It is noted that TrS* vanishes for the traceless algebra like that associated with SU{2)j, 
while gives dimension of the s.p. space for U{l)z and U{1)n. 
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In order to apply the generalized Wick's theorem to reduce the quantities in (jUj), we 
introduce an operator symbol [Q] [l6^, and express the TFD vacuum expectation value of 
an operator O as 

(0T|6g|0T) = (OT|g|OT)(OT|o[g]|OT) ; [Q] = --2-—. (59) 

Then, the basic matrix elements generalized by the TFD, which are building blocks for the 
quantum number projection at finite temperature, are provided by 

A, , ^ {QT\mlPl\QT) = [MK-\u = ~[{K'T'M'%,, (60a) 

B, . = {Ot\MAQ]\Ot) = [K-'NU = -[N''{K'T\u, (60b) 
C^u = {OrrnQHlOT) = [K^\u, (60c) 

where the suffices /i and u represent all the possible q.p. states. The TFD version of the 
generalized Wick's theorem |l6| is exemplified by 



{OT\miPiPlPi\OT) 


— A A 


A A 

■'^p.p-'^va 


ApfjAjyp^ 


(61a) 




— r A 


C A 

^ pp-'^ua 


-J- r" A 

~ ^pCT-'^Up) 


(61b) 


{0T\(3MQ]f^X\0T) 


— fl A 


C ppC ycr 


~l~ C p(jC up^ 


(61c) 


{Ot\MuWH\Ot) 


Bfiu C^pa 


B p^pC ua 


~l~ CpfjBi/pj 


(61d) 




Bpu Bpcr 


B ppByfj 


~l~ B p^Bup. 


(61e) 



For the SU{2)j x U{l)z x U{1)n projection, we have 

h^S = (^pfip + y^A, (62) 

where the quantity Qp = J2jep(3 + V^) (^n = J2jenU + V^)) represents a half- number of 
total s.p. levels in the model space. Hence, we obtain 

f27r 1 r2TT 



2 7 4-1 r'^'^ u) I"" r'^'^ 

' sm^-duj / sin^rf^ / rf0xj(^)(detK)i/2^ ^gg^ 
2 Jo Jo 



An 
and 



27r Jo 271 Jo 

2 7+1 f"^^ id f"^^ 
X ^ , / sin^-rfu; / sinOde / d(l)Xjiuj){detKy^^ {Ot\H[Q]\Ot)- 
Jo 2 Jo Jo 



(64) 
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Applying Eq. to (Or|if[Q] |0t), we can compute the Hamiltonian kernel in Eq. (j64j) . The 
SU{2)j X U{\)z X U{1)]^ projected thermal energy is calculated from Eq. by carrying 



out integrals in Eqs. ()63l64p . It is pointed out that the method o: 



integration, which has 



been demonstrated to~....^l.l\. the zero-temperature case g i, . avarlable also for 

the projected statistics. If we choose Euler angles (a, /?, 7) rather than the variables (cj, 6*, 0) 
for integral variables, the Gauss- Chebyshev quadrature formula will be useful for the a and 
7 integrals as well as the yjp and tp^i integrals, and the Gauss- Legendre quadrature formula 
for the j3 integral, since the function Xj{'^) is nothing but a sum of the Wigner D-functions. 

It should be noticed that the original q.p. occupation number is affected by the projection. 
Using the TFD transformation in (|T8|l . we get 

ala, = f, + {l- f,)f3lP, - ^JUl-f,)PA + ^U^-fM^i- (65) 

The first and third terms in the above expression contribute to the projected occupation 
number, when the projection operator is placed on the right of aj^a^. For P[j^z,n), the 
projected q.p. occupation number is calculated by 

^TT Jq 2tt Jq 

2J + 1 f^^ . ^Lo^ r . _ f^\^ 

X -— / sm —auj / smudO / am 



xxjH(detK)V2[/^ _ ^f^{l-f^){K~'N),f,]. (66) 

The present formalism is summarized by Eqs. (j46l58l59l6ip for general cases, with A^,^, 
B^v, Cfiu of Eq. (jUUI) . and K^^, = Z*^, M^,^ = iV*^, of Eq. (jHH), which are functions of the 
GBT coefficients {Ukfj.,Vk^,U^^,V,^^} and the q.p. occupation numbers The integral 
in Eq. ^ is exemplified by Eqs. (I6;^I64I66I1 for the SU{2)j x U{l)z x U{1)n projection. 
The resultant formulae are equivalent to those derived in Ref. jl^], despite the difference in 
appearance. It still deserves noting that the present formulae are straightforward extension 
of the zero-temperature formulae of Refs. 



D. Representation in specific bases 



For practical use, we here present explicit forms of the quantities relevant to the SU{2)j x 
U{l)z X U{1)n- projection, taking specific s.p. bases. 

1. Spherical bases 

Each s.p. level \Tnljm) is specified by r(= p,n), radial quantum number ra, orbital 
angular momentum /, total angular momentum j (= /± 1/2) and its projection to the 2;-axis 
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m. Here the s.p. state k is regarded as an abbreviation of {Tnljm), namely, \k) = \Tnljm). 
Analogously, \k') = \T'n'l' j'm') . For the projection P(j^z,n), we have Q = e"*"^^/?, yielding 
the Q-matrix in Eq. (j48j) of 

Qkk' = Q ml jm,T' n'l'j'm' = ^TT'^nn'^ll'^jj' ^^'^^ ^rn'm' i^'^) 

The K, L, M and N matrices are explicitly written as 



K — T* 



fill 



g{V^'e'^-D^U + U^'e-'^-D^'V)g g{V^' e'"^- DW* + U^'e''^^ D^'U*)g + f 

(68a) 

^ _ _ . g{V^'e'^-D^U + W'e~'^-D^'V)g g{V^'e''^- D^V* + We'^^ D^U*)g - gg\ 

~ ~ ^ g(We'^^D^U + V^e-'^-D''V)g -gg g(We'^- D^V* + V^e-'^^ D^'U*)g j ' 

/ /if 

(68b) 

where, e.g., 

[giU^e'^P^D^U + Vh^'^^^ D''V)g + g\u 
= E ^rr'5nn'5ii'5,,, g .{pi^e^''^ D^.^U,,, + V^^c-^^^ D^J/,, ^)g^ + gl5^^. (69) 

k,k' 

When we take only the particle number projection, all -D^,^ should be replaced by 5mm'- 
2. Goodman's bases 

Whole spherical basis states for a given set of quantum numbers [rnlj) can be separated 
into two sets {\Tnljm)} and {\Tnljm)}. The first set {\Tnljm)} consists of the states 
having (m — 1/2) equal to even integers, while the second set {\Tnljm) = T\Tnljm)} having 
(m — 1/2) equal to odd integers, where T is the time-reversal transformation. The phase 
convention is assumed to be 

\Tnljm) = T\Tnljm) = {—iy^"^^^\Tnlj — m), T^lmljm) = —\Tnljm). (70) 

Since the auxiliary Hamiltonian in Eq. (0} is invariant under the rotation about the x-axis 
by angle vr, i.e. = e"*'"'^'', the solution of the CHFB equation usually becomes eigenstates 
of R^. Therefore, a convenient choice is to employ the set of eigenstates of R^ for s.p. basis. 
From the relations as 

^^^Tnljm^x ~ ^ '^^rnlj-mi i^x) '^rnljmi^x) ~ ~^Tnljmi C^-'-) 

we find that the eigenvalues of R^ are ±i, 

Rj^Rl = -tcl rArI = +icI, (72) 
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where the corresponding eigenoperators are given by 

' _lY+m+l/2 



The operator i^j, is called signature operator, and the bases (kk) are called Goodman's 
bases j2^. It is sufficient to adopt the bases of Eq. (|7!?|l only for even (m — 1/2) to avoid 
repetition, which will symbolically be indicated by /c > 0. Using the relations in Eq. (|70|l . 
we can confirm that ct and c|; are related to each other by the time-reversal transformation. 



TclT 



TciT 



— c 



- (74) 

Accordingly matrix forms for expectation values of relevant quantities are also simplified. 
As an example, we here present the angular momentum components. 



(kk),(k'k') 



{I 



z)(kk),(k'k') 



{k\J,\k') 

{k\J.\k') 

{k\J,\k') 
{k\Jz\k') 



{Jy, 



{kk),{k'k') 



{k\Jy\k') 
{k\Jy\k') 



(75) 



Due to time-reversal properties, we have the following relations, 

{k\l\k') = -{k'\l\k), {k\jy\k') = {k'\jy\k), {k\l\k') = {k'\l\k) . 



(76) 



Unless the signature is spontaneously violated, the q.p. operators aj^ and aj^ belonging 
to the HFB eigenvalues E^^ and E/i, respectively, are classified in a similar manner. 



R^ajiKl = +iali and Ta]^T 



a 



A" 



(77) 



The conservation of the signature simplifies the structure of the GBT matrix as 



^k 

\4/ 



^C-A fa,\ 




u,-, 

V,-, Ul^ 



(78) 



\4/ 



Therefore, the TFD-extension of the GBT is represented as 

/ Cfc \ 



^{kk) 
- 

^{kk) 



H>0 



(fcfc) 



Cfc 
VCfc/ 



(mm) 



V4/ 



(79) 
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The transformation matrix W is given by 

(fcfc), ififi) It/- tt* _ / ' 

\ ^{kk),{flil) ^(kk),(tMp.) / 

I U,,g, VCf,g-, \ 

U-,,g, V-*^g, 



V(kk), {,,p.) 



" "^kfiyfi ''k^l^^^ ^ 

/ Vum Ut,g^. \ 
V-u,-g, 

-f/,V(7, Vk-,-g-, 



(80a) 



(80b) 



We now consider the SU{2)j x U{l)z x ?7(l)Ar projection. The operator Q = e ^'^^R 
gives the transformation rules for the s.p. operators 

QckQ^ = e*^- J2 f^fcVcfc' + Glf:,c^}j , Qc^gt = e^^. ^ ^^^,,0,, + F^^.c,-,] , (81) 

k'>0 ^ ^ fc'>0 ^ ^ 



where the coefficients F and G are defined by 

1 



Fkk' — Fkk' — ^TT'^nn'^ll'^jj' 7^ 



m' ,m ' —m',—m 



lV-i/2/'r)i I n-? 

/ V m' ,—m ' —m'.r. 



(82a) 



j^_]_y+m'+l/2 



m' ,m —m' ,—m ' V / V m',— m ~m' ,m) 



(82b) 



The transformation ()8ip is converted to the one among the TFD q.p. operators via 
Eqs. dZMD, 



Q I ^\^^^ ] = 5Z ^ ^^^^^^ ^''^^ ^^^^'^^ 



(83) 
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where the submatrices K, L, M and N are provided by 



g-^{We'^^GU + V^e-'^^GV)-^,g, gp{We'^-FU + e'^^^ FV)-^,g, 
g^{ytr^ivrQu + U^'e-'^-GV)^^g^ g^{V^'e''^- FU + U''e-'^^ FV)^pgp 
\ g-^{y''e'^^F*U + U''e-'^^F*V)j,^g^ g^{V''e'^^G*U + C/*'e-'^-G'*\/)^p^^ 

g^{We'^-G*V* + V^e-'^^G*U*)^,g, g^{U^ e'^^ F*V* + e-'^^ F*V*)^^g^ \ 
gpiWe'^^FV* + V^e-'^^FU*)f,,g, g^iWe'^^GV* + V^e-'^-GU*)-^^g^ 
g^iV^'e'^-FV* + W'e-'^- FU*)^^,g^, g^,{V^'e'^-GV* + U^'e-^^-GU*)^pgp 

g.{y''e'^^G*V* + U^'e-''^^G*U*)-^^g^, g-^{V^'e'^^ F*V* + W'e-'^^ F*U*)ppgp ) 



M, 



gl 

^2 

V ^^y 

m 



(84a) 



( gf,{V''e"^-GU + U''e-''P-GV)^^g^ g^^iV^'e'^^ FU + We-''^^ FV)^pgp 
gpiV^'e'^-F*U + U^'e''^- F*V)p^g^ g^{V^'e'-^-G*U + U^'e-'^^G*V)ppgp 
g^{We'^-F*U + V^e-''^^F*V)^,g, g^{We'^-G*U + V^e-'^-G*V)^pgp 
V gfi(We"^-GU + V^e-"^^GV)fi,g, gi,{We"^^FU + V''e-"^^FV)ppgp 

g^{V^'e'^-FV* + W'e-'^- FU*)^^g^ g^{V^'e'^-GV* + U^' e-^^-GU*) ^j^g^ \ 
gj,{Ve'^^G*V* + U^'e-''^^G*U*)-^^g^ gp{V^'e'^- F*V* + W'e-'^^ F*U*)^pgp 
g^(U'^e''^-^G*V* + V'^e-''^^G*U*)^^g^ g^{U''e''^^ F*V* + V^ie'''^^ F*U*)f,pgp 
gpiWe'^^^FV* + V'^e-'^^FU*)~^,g, gpiWe'^^^GV* + V^e-"^^GU*)i,pgp J 

/ g^g^ \ 
g^gp 



- 5 



fJ,U 



g/^gti 

V gugn 










/ 



(84b) 



Note that, for the case of particle number projection only, all the matrices F should be 
replace by an identity, and G by zero. 

E. Formulae in BCS approximation 

We continue to work with the Goodman's bases. The HFB approximation is reduced to 
the BCS approximation by the replacements 



u-> 



fe/i 



(85) 
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We here assume Uk and Vk to be real. Thereby the BCS transformation is expressed as 



/ Cfc \ 



V4/ 



/ Mfc Vk\ f ak\ 

Uk -Vk 

Vk Uk 

\-Vk UkJ \alJ 



a. 



(86) 



Accordingly, the quantum number projection is somewhat simplified in the BCS approxi- 
mation. By using the notations 



^kk' = UkUkie'-f-' + VkVkie '"^^ = {ukUk' + VkVk') cos(/?^ + i{ukUk' - VkVk') sint/?^, 

rikk' = UkVk'e"^^ - VkUk'e"'"^^ = {ukVk' - VkUk') cos if r + i{ukVk' + VkUk') sin ipr, (87) 

the matrices K, Z, M and N for the SU{2)j x U{l)z x C/(l)Ar projection are obtained by 

-^(fefc), (fc'fc') = L\kk), (k'k') 

( 9k^kk'F^k'9k' 9kikk'G*j^^,gk, —gkilkk'G*j^^,gk' gkVkk'F^k'dk' \ 



dkikk'Gkk'gk' dkikk'Fkk'dk' 

9kV*kk'Gkk'9k' 9kVkk'Fkk'9k' 

V -9kVkk'Fkk'9k' -g-kilkk'G*kk,9k' 

fgl \ 



-9kVkk>Fkk'9k' ■~9k''lkk'Gkk'9k' 
9kikk'Fkk'9k' 9kCkk'Gkk'9k' 
~9kikk'G*kk,9k' 9kikk'^kk'9k' / 



kk' 



(?! 

^2 
V Q -glJ 



(88a) 



{kk), {k'k') 



N, 



{kk), {k'k') 



( -9kVkk'Gkk'9k' 
-9kVkk'Fkk'9k' - 
9k^kk'F^f^,gk' 

\ —9k^kk'Gkk'9k' 

/ 


gk9k 

V gkgk 



>kk' 



gkVlk'Fkk'gk' 

-9kVkk'Gl,^,gk 
gkikk'G*j^j:,gk, 

dkikk'Fkk'dk' 

QkQk \ 

gm 



y 



9k^kk'^kk'9k' 
-9kikk'Glj:,gk' 

-9k1lkk'Gl,j^,gk' 
—9kVkk'Fkk'9k' 



9kikk'Gkk'9k' \ 
9kikk'Fkk'9k' 
9kVkk'F^i^,gki 

-9kVkk'Gkk'9k' ) 



(88b) 
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For the case of the number projection only, putting F = 1 and G = in the above expression 
and employing ^kk = cosy^r + i{ul - vl) sin^j^, 77^^ = -r]l^ = 2iukVk sin ip^, we get 



K(kk),(k'k') - ^kk),(k'k') 



^kk' 



f89a) 



M(^kk),(k'k') — ^{kk),{k'k') 



Skk' 



I Mkk + gl 9k9-kVkk \ 

glikk + gl -9k9kVkk 

-QkQkVkk gUkk + gl 

V gmvkk gi^i^ + gij 

( -gkgp]kk gkgkiilk - 1) \ 

grnvkk gmi^kk - 1) 

gkgkiikk - 1) gkg-k^kk 

V gmiikk - 1) -grnvkk o / 

(89b) 

Eq. M^ yields 

{deiKf/^ = n [fk{l - f-k) + (1 - + (1 - fk){l - fk)ikk + fkhCkk] ■ (90) 

fc>0 



Eqs. (I89l9(jj) give the formulae equivalent to those in Ref. 
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IV. APPROXIMATION OF ENTROPY 

Once the original Bogoliubov transformation and the q.p. energy in Eq. ^ are 
determined, the projection method discussed in Sec. lllll is straightforwardly applicable. This 
indicates the variation-before-projection (VBP) calculations. However, there still remains a 
serious problem in formulating the variation-after-projection (VAP) calculations [l^ . 

Between the exact free energy p'^^^'^^^ which is generated by the exact Boltzmann-Gibbs 
operator exp(— i^/T), and an approximate one Fp, which is generated by P exp{—HQ/T)P 
in our case, the Peierls inequality holds [l(^ . 



-rexact 



-Tln[Tr(e"^/^P)] < Fp = Ep - TSp. 



(91) 



This inequality vindicates the variational calculations minimizing Fp. The approximate 
energy Ep is expressed in terms of the TFD vacuum expectation value by 



Ep = Ty{wpH) 
The entropy is given by 

Sp = — Tr(wplnti)p) = — Tr 



Tt{woHP) _ {Ot\hP\Ot) 
Tr{woP) (Or|P|Or) 



(92) 



iJo/Tp Y 
Tr(e-^o/Tp) - \^ Tr(e-^o/^P) J 



-In 



(93) 
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where a relation Pln{PAP) = Pln{AP) for any operator A is applied. However, no further 
reduction is allowed because [P, Hq] ^ 0. Since the projection operator remains in the 
logarithmic function, it is prohibitively difficult to deal with the entropy in Eq. ()93|) without 
further approximation. 

In order to make the VAP calculations practical, we here discuss an additional approxi- 
mation on the entropy. The VAP scheme is based on the fact that an approximate entropy 
does not exceed the exact one, or the Peierls inequality (PT|). Therefore, any approximation 
method adopted in the VAP scheme should not conflict with those inequalities. 

We now introduce a projected space V and denote its complementary space by Q. The 
operator Pe~^°^^P is hermitian and non-negative. We here denote eigenstates of Pe~^°^^P 
by \K) and its eigenvalue by ujk{> 0), 

Pe-^o/Tp \k) = \K)ujk- (94) 

Any states in the Q are eigenstates of Pe~^'^^^P with the null eigenvalue. The V space is 
also spanned by eigenstates of Pe~^°^^P. For \K) eV, P\K) = \K) leads to 

UK = (ir|Pe-^«/^P|i^) = (ir|e-^«/^|ir) . (95) 

Therefore, an essential part of the entropy is expressed as 

Tr[Pe-^»/^Pln(Pe-^»/^P)] = cu,^ Iuuk = ^ (ir|e"^°/^|ir) ln{K\e-''°/^\K) . (96) 

Observing that the functional form of the above expression is the convex function F{x) = 
X Inx with X = {K\e~'^°^'^\K), we apply the inequality ()B2jl in Appendix B to derive a 
useful inequality as follows, 

Tr[Pe-^«/^Pln(Pe-^«/^P)] < (ir|e-^«/^ In e-^°/^|ir) 

Kev 

i Y: {K\Pe-^"!^H,\K) = -^T.{e-^^^I^H,P). 

KeV+Q 



(97) 



Thus, we obtain an inequality between two approximate entropies Sp and Sp, which is 
deflned below. 



Sp = -Tr 



Pe-^o/Tp / pe-Ho/rp 
Tre--f^o/^P \ Tre--^o/^P 



>lI^^(^:^^ + lnTr(e--/-P).5^. (98) 
- T Tr(e-^o/rp) K ) p \ J 
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The approximate free energy defined by Sp satisfies the Peierls inequahty, 



F^ = Ep- TS'p >Fp = Ep- TSp > F'^"'^"*. (99) 

It is noted that this approximate entropy Sp can also be obtained by deahng with the 
logarithmic function in Eq. ()93|) by Pln(e~^"/"^P) ~ —PHqP/T, as if Hq and P were 
commutable, though such an ansatz is not adopted in the present argument. An crucial 
point is that this approximation preserves the Peierls inequality (PU)). which gives a ground 
for the variational calculations. Therefore, Fp will be available for the YAP calculation, 
since the smaller Fp gives the better approximation to The TFD expression of Fp is 

given by 

= ^°"'^f"£;^f'°"^ - Tln(0.|P|0.) +tJ: hi - U). (100) 

In Fp, all the quantities can be calculated in the formalism presented in Sec. IIIII For the 
SU{2)j X U{l)z X U{1)n projection, (Ot|F|Ot) and {'^t\HP\Qt) are given by Eq. ^ and 
Eq. ()64|) . respectively. Moreover, (OT|ffo-P|OT) is immediately obtained from Eq. (jUBj) . Since 
these elements are expressed as functions of the GBT coefficients f/^^, V^^, f/^^, V^^ and 
= sin^t?^ {i.e. coefficients of the TFD-extended GBT), the variation of Fp with respect 
to these variables produces a closed set of independent VAP equations. 



V. SUMMARY AND DISCUSSIONS 

Although the mean-field theories play an important role in nuclear physics at finite tem- 
perature as well as at zero temperature, the mean-field solutions often violate conservation 
laws. In order to restore desired quantum numbers for thermally equilibrated many-body 
states of an isolated finite system like a nucleus, we have applied the thermo field dynamics 
(TFD) to construct a transparent and practical formalism of the quantum-number-projected 
statistics. In the TFD, an ensemble average of an observable is expressed by a TFD vac- 
uum expectation value, whereas the single-particle operator space is doubled by introducing 
tilded operators. Making use of this method, we have derived formulae for the projected 
statistics at finite temperature. As an example, we have shown explicit representations for 
the projection of angular momentum as well as that of particle numbers. A significant ad- 
vantage of this formalism is to keep a complete parallel with the zero-temperature case, so 
that the projection method and the corresponding computer code, which have been demon- 
strated to be successful, could be directly translated to the projection at finite temperature, 
apart from the variation. 

Our formalism presented in Sec. IIIII is, in principle, applicable to both the variation- 
before-projection (VHP) and the variation-after-projection (VAP) schemes. In the mean- 
field calculations without the projections, sharp phase transitions appear; for instance, there 
could be a discontinuity in heat capacity at the critical temperature. This is not realistic 
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for finite systems, in which quantum fluctuations of the fields often wash out the distinct 
signature of transitions. While such unrealistic signature remains by the VBP calculations, 
the YAP scheme is expected to smooth it out, including a significant portion of the fluctua- 
tions. However, the VAP calculations minimizing the projected free energy are impractical 
without further approximation, since the entropy does not take a simple form because of 
the non-commut ability of P with Hq. Thus, in Sec. IIVI we have discussed an additional 
approximation for the entropy, which is expected to resolve this problem. It should be em- 
phasized that the Peierls inequality, which is an important requisite justifying the variation, 
is proven to be satisfied in this approximation. 

We have adopted the HFB or constrained HFB (CHFB) self-consistent solution at finite 
temperature for the basic quasiparticle picture, and have developed the TFD formalism on 
it. The CHFB approximation is the most effective mean- field theory in describing nuclear 
structure at low temperature or along the yrast. It is commented here that, whereas the 
constraints might not look important when the projection is implemented, it can still be 
useful. In the VBP calculation, it could be essentially important to obtain a mean-field 
solution well connected to the true many-body solution. Also in the VAP calculations, 
the CHFB solution may be used to obtain good initial configurations. Furthermore, if 
there remain some quantum numbers unprojected, Cj (j = 1,2, ■■■ ^Nc)^ they could be 
handled in the CHFB scheme in the projected statistics for other quantum numbers, by 
replacing the Hamiltonian H by the auxiliary one H' = H - Y.'.^-^ XjCj in the free energy 
in Eqs. (Wmm- 

The present formalism will be useful in investigating effects of quantum fluctuations 
connected to the conservation laws on thermal properties of nuclei. A topical example 
is the pairing correlations at finite temperature. Although the superfluid-to-normal phase 
transition has been predicted to occur at relatively low temperature (0.5 ^ < 1 MeV) [3, 
24, I, it is not yet well established in nuclei. Based on a precise measurement of nuclear 



level den.Me., the .upe.fluM-to-„o.mal t.an.Mon nuclei ha. been di.cu..ed recently Q. 
While there is no sharp discontinuity as in infinite systems, the S'-shape behavior in the 
graph of heat capacity vs. temperature, C(T), has been argued to be a signature of the 



transition 



26, 



27l |. For proper understanding of the phase transition and its relation to 
the S'-shape, it is important to investigate effect of the quantum fluctuations in connection 
with the conservation laws, although the other quantum fluctuations cannot be neglected 
in quantitative description j2S| , particularly around the critical temperature. This situation 
holds also for the deformed-to-spherical shape phase transition. 

Several theoretical frameworks have been invented to take into account the quantum 
fluctuations not restricted to those connected to the conservation laws; for instance, the 
static-path approximation (SPA) 0, H and the shell-model Monte Carlo (SMMC) 
approach ^jj. Whereas we have mainly focused our discussion on the projections in the 
mean-field approximations, it will be straightforward to incorporate the projections into the 
above extensive methods, in which the wave functions are represented by a superposition of 
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those in the mean-field theories. In practice, the SPA calculations with the particle-number 
and angular momentum projections have been carried out 0,E3|' based on the conventional 
thermal formalism. With its simplicity the present TFD formulation of projections will be 
useful also in this course. T he p roiections combined with the SMMC have been restricted 
to relatively simple cases l3l| . in which [P, ^o] = holds. The present formulation may 
help applying more general projections to the SMMC implementation; e.g. the angular 
momentum projection, and the particle-number projection in the pairing decomposition. 
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APPENDIX A: QUASIPARTICLE REPRESENTATION OF HAMILTONIAN 

We here present the quasiparticle representation of the Hamiltonian H' given in Eq. 
The Hamiltonian is assumed to consist of up to the two-body interactions as in Eq. Q. 
Applying the Bogoliubov transformation in Eq. (^, we rewrite H' in terms of quasiparticle 
operators. 



(vac.|i5f'|vac.) — \Z — X^N — u^ot 




(Alb) 



(Ala) 



(Ale) 




(Aid) 




(Ale) 




(Alf) 
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where 

{Hii)^, = [f/^(e + r(°))t/-yt(e + r(°))*y + f/tA(°V-y^A(°)*f/]^^, (A2a) 
{H2o)^.u = [f/t(^ + r(°))\/*-\/^(e + r(°))*f/* + f/^A(°)[/*- V^A(°)*\/*]^^, (A2b) 

ijkl 

+U:^V,:Uk.V,p - UiV,;U,^V,,] , (A2c) 



6 

ijkl 



+U;,UlU^„V,l + U:^V,,V:,V,l] , (A2d) 
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ijkl 



-UlU*,V:^V,l - U:^U;^V:,m] . (A2e) 

Definition of ^ is tlie same as given in Eq. (II2D. r(°) and A^o) are defined in an analogous 
manner to Eq. p3|) . with the single-particle density p^^-* and the pair tensor k^^^ defined by 

Plf = = (vac.lcjQlvac.) = iV*V'% , (A3a) 



= = (vac.|c,-c,|vac.) = (\/*f/*Oij • (A3b) 



APPENDIX B: PROOF OF INEQUALITY USED IN SECTION HVl 

Let F{x) be a real function, which is convex downwards (or concave upwards) in a simply 
connected domain of the variable x (> 0), and xq belongs to the same domain. This leads 
to the inequality 

Fix) > F{xo) + F'{xo){x - xo) . (Bl) 

As far as F"{x) > 0, the the left-hand side (Ihs) and the right-hand side (rhs) of Eq. ()B1|) 
become equal only at x = xq. For a hermitian and non- negative operator O and a given 
physical state the following relation turns out, 

{md)\z) > Fimz)) . (B2) 

This inequality is proven as follows. 

First, we introduce a complete set of the eigenstates of the operator O, denoted by {[A;)}. 
Each eigenstate \k) satisfies 

d\k) = \k)ujk. (B3) 

Then, the state \i) is expanded by 

\t) = J2\k)uk^, (B4) 

k 
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which is a unitary transformation. Using Eqs. ()B1IB4|1 . we obtain 

k 

> J2 [F{{z\d\z)) + F\{z\d\z))iuk - mz))] = F{{z\d\z)), (B5) 

k 

by applying Eq. ()B1|) with setting x = Uk and xq = {i\0\i), and using J2k l^fciP = 1 and 
J2k I'^kil'^ ^k = It is now obvious that the Ihs and the rhs of Eq. ()B2|) are equal only 

if 1^) is an eigenstate of O. 

It has been verified that, if F'(x) < and F"(x) > 0, we have 

5^(.|f(6)K)>5^f((.|6k)). (B6) 

i i 

We have here derived the inequality ()B2|) in more general manner, by lifting the sum over i 
and the condition F'{x) < 0, besides that truncation scheme is also discussed in Ref. \v^ . 
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